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HIGHLIGHTS 


►  We  provided  a  3D  numerical  model  of  a  TEG  module  used  in  a  waste  heat  recovery  system. 

►  The  effects  of  temperature  difference  and  waste  gas  heat  transfer  coefficients  were  investigated. 

►  The  power  density  for  a  TEG  module  is  the  objective  function  to  be  maximized. 

►  The  numerical  data  are  in  good  agreement  (within  8%)  with  the  experimental  data. 

►  The  optimum  TEG  module  spacing  and  spreader  thickness  are  strongly  dependent  on  the  waste  gas  heat  transfer  coefficient. 
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When  thermoelectric  generator  (TEG)  modules  are  attached  to  a  rectangular  chimney  plate  for  venting 
hot  flue  gases,  the  power  generated  per  unit  surface  area  (power  density)  is  strongly  dependent  on  the 
TEG  module  spacing.  The  thermoelectric  module  consists  of  a  hot  plate,  a  spreader,  a  thermoelectric 
generator  and  a  cold  plate  based  on  water  cooling.  In  this  study,  the  optimization  of  TEG  module  spacing 
and  its  spreader  thickness  as  used  in  a  waste  heat  recovery  system  is  investigated  and  solved  numerically 
using  the  finite  difference  method  along  with  a  simplified  conjugate-gradient  method.  The  power 
density  for  a  thermoelectric  module  is  the  objective  function  to  be  maximized.  A  search  for  the  optimum 
module  spacing  (S)  and  spreader  thickness  (Hsp),  ranging  from  40  mm  <  S  <  300  mm  and 
1  mm  <  Hsp  <  30  mm,  respectively,  is  performed.  The  effects  of  different  operating  conditions,  including 
the  temperature  difference  between  the  waste  gas  and  the  cooling  water  (AT  =  200—800  K),  and 
effective  waste  gas  heat  transfer  coefficients  (fy,  =  20-80  W/m2  I<)  are  discussed  in  detail.  The  predicted 
numerical  data  for  the  power  vs.  current  (P— /)  curve  are  in  good  agreement  (within  8%)  with  the 
experimental  data. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Thermoelectric  generators  present  potential  applications  in  the 
conversion  of  low  level  thermal  energy  into  electrical  power. 
Especially  in  the  case  of  waste  heat  recovery,  it  is  unnecessary  to 
consider  the  cost  of  the  thermal  energy  input,  and  there  are  addi¬ 
tional  advantages,  such  as  energy  saving  and  emission  reduction,  so 
the  low  efficiency  problem  is  no  longer  the  most  important  issue 
that  we  have  to  take  into  account  [1  ].  In  general,  a  TEG  consists  of 
a  number  of  semiconductor  pairs  that  are  connected  electrically  in 
a  series  and  thermally  in  parallel,  and  each  pair  includes  a  p-type 
and  an  n-type  element.  Although  in  theory,  a  single  piece  of 
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semiconductor  material  could  work,  a  series  connection  is  used  to 
meet  the  high  voltage  potential  requirements.  p-Type  and  n-type 
elements  are  alternated  to  assure  that  the  carriers  transport  in  the 
same  direction. 

Most  of  the  reported  performance  data  of  TEGs  is  analyzed  using 
conventional  non-equilibrium  thermodynamics  [2,3].  With  refer¬ 
ence  to  the  structure  of  a  TEG,  a  significant  increase  in  the  power 
output  from  a  module  can  be  achieved  by  modifying  the  geometry 
of  the  thermoelectric  elements  [4,5].  In  subsequent  research,  Rowe 
[6,7]  provided  efficiency  in  a  couple  of  solar-powered  TEGs  and 
reviewed  applications  of  nuclear-powered  TEGs  in  space.  Chen 
et  al.  [8]  used  an  irreversible  model  to  study  the  performance  of 
a  TEG  with  external  and  internal  irreversibility.  The  optimal  range 
of  the  parameters  for  the  device-design  was  determined,  and 
the  problems  relative  to  the  maximum  power  output  and 
maximum  efficiency  were  discussed.  Xuan  et  al.  [9]  employed 
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Nomenclature 

Greek  symbols 

a 

Seebeck  coefficient  (V/K) 

A 

TEG  module  area  (mm2)  A  =  S  x  S 

£ 

descend  direction  coefficient 

Dh 

hydraulic  diameter  (m) 

£ 

flue  gas  emissivity 

E 

electric  field  (V/m) 

V 

TEG  conversion  efficiency  (%)  p  =  PmaxlQh 

h 

effective  heat  transfer  coefficient  (W/m2  K) 

0 

general  dependent 

H 

height  (mm) 

P 

electrical  resistivity  (Q  m) 

I 

electric  current  (A) 

a 

Stefan  Boltzmann  constant  (W/m2  K4) 

j 

electric  current  density  (A/m2) 

? 

search  directions 

Jobj 

objective  function 

k 

thermal  conductivity  (W/m  K) 

Subscripts 

L 

length  (mm) 

c 

cold  side 

n 

pairs  of  p-type  and  n-type  semiconductor 

cer 

ceramic  substrate 

P 

power  output  (W) 

cond 

conductive  copper 

Pmax 

maximum  power  output  (W) 

conv 

convection 

P  maxi  A 

maximum  power  density  (W/m2) 

f 

fluid 

d 

heat  transfer  rate  (W) 

gas 

flue  gas 

R 

residual  sum 

h 

hot  side 

S 

TEG  module  spacing  (mm) 

hp 

hot  plate 

T 

temperature  (I<) 

opt 

optimal  design 

AT 

temperature  difference  between  the  waste  gas  and  the 

P-N 

p-  and  n-type  thermoelectric  semiconductor 

cooling  water  (I<) 

rad 

radiation 

V 

electric  potential  (V) 

sp 

spreader 

V 

velocity  (m/s) 

TEG 

thermoelectric  generator 

Voc 

open  circuit  voltage  (V) 

water 

cooling  water 

*i 

design  variable  vector 

x,y,  z 

Cartesian  coordinates 

a  phenomenological  model  to  study  the  effects  of  internal  and 
external  interface  layers  on  thermoelectric  module  performance. 
Sets  of  general  performance  formulas  have  also  been  derived. 
Pramanick  and  Das  [10]  performed  a  study  on  the  structural  design 
of  a  thermoelectric  module.  A  model  for  a  cascaded  TEG  was 
developed,  which  was  based  on  finite  time  thermodynamics. 

Recently,  there  has  been  a  growing  interest  in  TEGs  using 
various  heat  sources  in  such  applications  as  combustion  of  waste, 
geothermal  energy,  power  plants,  and  other  industrial  heat¬ 
generating  processes  [11,12].  Khattab  and  Shenawy  [13]  studied 
the  possibility  of  using  a  solar  TEG  to  drive  a  small  thermoelectric 
cooler  (TEC).  Crane  and  Jackson  [14]  studied  numerical  heat 
exchanger  models  integrated  with  models  for  Bi2Te3  thermo¬ 
electric  modules  which  were  validated  against  experimental  data 
from  previous  cross-flow  heat  exchanger  studies  as  well  as 
experiments  using  thermoelectric  modules  between  counter-flow 
hot  water  and  cooling  air  flow  channels.  Dai  et  al.  [15]  provided 
a  new  type  of  TEG  system  based  on  liquid  metal  which  served  to 
harvest  and  transport  waste  heat.  The  experimental  results  for  the 
TEG  system  were  discussed  and  a  calculated  efficiency  of  2%  in  the 
whole  TEG  system  was  obtained.  Champier  et  al.  [16]  reported 
a  study  conducted  in  order  to  investigate  the  feasibility  of  using 
a  TEG  in  an  improved  biomass  fired  stove.  The  maximum  power 
reached  by  each  module  varied  between  1.7  W  and  2.3  W, 
resulting  in  a  temperature  difference  between  the  two  sides  of 
160  °C. 

Furthermore,  in  the  application  of  TEGs  for  waste  heat  recovery 
power  generation,  there  have  been  many  conceptual  designs  for 
a  power  conversion  system  that  are  potentially  capable  of  being 
applied  in  this  area.  Hsiao  et  al.  [17]  constructed  a  mathematic 
model  to  predict  the  performance  of  a  TEG  module  attached  to 
a  waste  heat  recovery  system.  The  results  showed  that  the  TEG 
module  presented  better  performance  on  an  exhaust  pipe  as 
opposed  to  a  radiator.  Thacher  et  al.  [18]  investigated  the  feasibility 
of  waste  heat  recovery  from  exhaust  in  a  light  truck  by  connecting 


a  series  of  16  TEG  modules,  which  showed  good  performance  at 
high  speeds.  Niu  et  al.  [19]  constructed  an  experimental  TEG  unit, 
which  was  used  to  examine  the  influences  of  the  main  operating 
conditions,  the  hot  and  cold  fluid  inlet  temperatures,  flow  rates 
and  the  load  resistance  on  the  power  output  and  conversion  effi¬ 
ciency.  Karri  et  al.  [20]  studied  the  power  and  fuel  savings 
produced  by  TEGs  placed  in  the  exhaust  stream  of  a  sports  utility 
vehicle  (SUV).  The  optimized  quantum- well  (QW)  based  TEG  stack 
generated  about  53—5.8  kW,  resulting  in  a  fuel  savings  of 
about  3%. 

According  to  the  literature  survey  [21  ],  it  is  recognized  that  most 
of  the  previous  theoretical  models  have  been  limited  to  one 
dimensional  problems.  However,  in  these  existing  models,  the  P—N 
element  pair  is  simply  treated  as  a  single  bulk  material  so  that  the 
difference  in  thermal  behavior  between  the  two  semiconductor 
elements  was  not  possible  to  evaluate.  Only  a  few  studies  have  been 
conducted  for  three-dimensional  TEGs.  For  example,  Chen  et  al. 
[22]  proposed  a  three-dimensional  TEG  model  and  implemented  it 
in  a  computational  fluid  dynamics  simulation  environment 
(FLUENT)  and  its  user-defined  functions  (UDFs).  Cheng  et  al.  [23] 
provided  a  three-dimensional  numerical  model  to  predict  tran¬ 
sient  thermal  behavior  of  TECs.  It  was  observed  that  both  the 
temperatures  of  the  hot  and  the  cold  ends  increased  with  the 
cooling  load,  and  the  value  of  the  COP  (coefficient  of  performance) 
linearly  increased  with  the  cooling  load.  Recently,  Sidek  et  al.  [24] 
provided  a  numerical  study  on  temperature  uniformity  of 
thermal  spreader  integration  with  microcombustor  for  TEG.  The 
results  show  the  effect  of  thermal  conductivity,  specific  heat 
capacity  and  thickness,  and  the  mutual  relation  between  the 
magnitude  of  heat  source  and  thermal  spreader  temperature 
profile.  Jang  et  al.  [25]  used  a  numerical  optimization  technique  in 
a  geometrical  optimization  for  spreaders  to  obtain  efficient  power 
output  balancing  with  reasonable  material  costs.  The  results 
showed  the  optimum  spreader  length  and  thickness  with  the  fixed 
module  spacing  for  Pareto-optimal  configurations. 
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In  the  present  study,  the  optimization  of  TEG  module  spacing 
and  its  spreader  thickness  is  investigated  and  solved  numerically 
using  the  finite  difference  method  along  with  a  simplified 
conjugate-gradient  method  (SCGM)  [26].  To  reach  optimization 
goals,  the  power  density  for  the  thermoelectric  module  is  the 
objective  function  to  be  maximized.  The  effects  of  different  oper¬ 
ating  conditions,  including  the  temperature  difference  between  the 
waste  gas  and  the  cooling  water  (AT  =  200-800  K),  and  effective 
waste  gas  heat  transfer  coefficients  (hh  =  20-80  W/m2  K)  are 
investigated  and  discussed  in  detail.  In  addition,  experiments  are 
conducted  to  verify  the  physical  model. 

2.  Mathematical  analysis  for  the  thermoelectric  module 

The  schematic  diagram  of  a  waste  heat  recovery  system  with 
thermoelectric  modules  is  shown  in  Fig.  1.  The  thermoelectric 
module  consists  of  a  hot  plate,  a  spreader,  a  TEG  and  a  cold  plate 
based  on  water  cooling.  The  dimensions  of  the  assembled  square 
TEG  unit  are  30  mm  x  30  mm  x  3.42  mm.  Its  material  is  B^Tes  and 
the  detailed  geometrical  parameters  and  thermal-electrical 


properties  of  the  TEG  module  are  listed  in  Table  1.  A  TEG  is  the 
major  component  of  the  thermoelectric  module,  in  which  31  pairs 
of  p-type  and  n-type  semiconductor  legs  composing  the  generator 
are  connected  thermally  in  parallel  between  the  hot  waste  gas  and 
cold  plate  and  electrically  in  a  series  to  power  the  load  circuit.  Since 
we  use  the  thermoelectric  models  on  the  basis  of  an  equivalent 
current  method,  which  means  there  is  an  equal  electric  current  and 
an  equal  cross  area  of  thermoelectric  elements,  only  one  pair  is 
considered  here.  Table  2  shows  the  equivalent  thermal-electrical 
properties  (thermal  conductivity  kp-u,  electrical  resistivity  pp-u 
and  Seebeck  coefficient  ap-u)  when  a  TEG  module  with  n  pairs  of 
p-type  and  n-type  semiconductors  is  simplified  to  one  pair,  where 
n  =  31  in  the  present  study. 

2.1.  Governing  equations 

Assuming  constant  thermal  and  electrical  properties  of  the 
thermoelectric  material,  the  governing  equation  of  the  thermal 
field  (T)  and  electrical  potential  (V)  for  the  p-type  and  n-type 
element  can  be  expressed  as 


Fig.  1.  Schematic  diagram  of  the  thermoelectric  module  with  a  spreader  for  waste  heat  recovery. 
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Table  1 

Geometrical  parameters  and  thermal— electrical  properties  of  the  TEG  module. 
Geometrical  and  thermal— electrical  properties  Symbol  Value 
1.  Hot  plate 


Thickness  (mm) 

Hhp 

30 

Thermal  conductivity  (W/m  K) 

khp 

43  (Steel) 

2.  Spreader 

Thermal  conductivity  (W/m  K) 

ksp 

155  (Aluminum) 

3.  Ceramic  substrate 

Thickness  (mm) 

Hcer 

0.75 

Thermal  conductance  (W/m  I<) 

her 

49.2 

4.  Conductive  copper 

Thickness  (mm) 

Hcond 

0.12 

Thermal  conductance  (W/m  I<) 

kcond 

403 

Seebeck  coefficient  (V/K) 

&cond 

1.4  x  10"5 

Resistivity  (Q  m) 

Pcond 

1.44  x  10"8 

5.  P—N  thermoelectric  elements  (for  one  pair) 

Thickness  (mm) 

Hp-N 

1.68 

Thermal  conductance  (W/m  I<) 

kp-N 

1.6 

Seebeck  coefficient  (V/K) 

ap-N 

±2.0  x  10~4 

Electrical  resistivity  (Q  m) 

Pp-n 

9.0  x  lO"6 

kP_N 


d2T 

ay2 


_^2 

+  Pp~nJ  =  0 


(1) 


where  hh  is  the  effective  heat  transfer  coefficient  combing  the 
radiation  ( hrad )  and  convection  (hCOnv)  heat  transfer  coefficients, 
which  may  be  estimated  as  follows  [27]  : 

hrad  =  0-em(Ts  +  Tgas )  •  (t2  +  T2a^j  (5a) 

Nu  =  0.023  Re0  8Pr1/3  =  ^con^  Dh  Dittus-Boelter  equation 

(5b) 

where  a  is  the  Stefan  Boltzmann  constant,  and  s  is  the  flue  gas 
emissivity.  The  Reynolds  number  is  Re  =  vDh/v  and  Pr  is  the  Prandtl 
number.  Dh  and  kf  are  represent  the  hydraulic  diameter  and  thermal 
conductivity  of  the  fluid,  respectively.  It  can  be  shown  that  for  Tgas 
ranging  from  500  to  1100  K  and  flue  gas  velocity  ranging  from  1  to 
5  m/s,  the  effective  heat  transfer  coefficient  hh  is  ranging  from  20  to 
80  W/m2  K. 

(2)  At  the  interface  between  the  hot  plate  and  spreader,  there  is 
a  continuity  of  temperature  and  heat  flux. 


Tsp  =  Thp  and 


(6) 


(3)  At  the  interface  between  spreader  and  ceramic  substrate. 


E 


~  aP-N 


dTy  dT -r 

3x '  +dyJ  + 


Trer  —  T< 


sp 


dT 

ray 


<ind  —  hcer^T,  —  -ksp— 


dT 

’ey 


(7) 


-  Pp-nJ 

(2) 

where  kp- jy,  pp-jv  and  ap _jv  represent  the  equivalent  thermal 
conductivity,  electrical  resistivity  and  Seebeck  coefficient  of  TEG 
module  with  31  _pairs  of  p-type  and  n-type  semiconductors, 
respectively.  The  E  and  J  denote  the  electrical  field  and  current 
density,  respectively. 

The  energy  heat  conductions  for  the  solid  spreader,  ceramic 
substrate  and  conductive  copper  are 


(4)  At  interface  between  ceramic  substrate  and  conductive 
copper. 

t  t  a  i  dT  .  dT  (  . 

her  —  hond  ana  —  <cer0y  —  ~  <cond^y  W 

(5)  At  the  interface  between  conductive  copper  and  p-type  and 
n-type  thermoelectric  elements, 


a2r  a2r  a2r  . 

— Y  H - - n  =  0 

Ox2  Oy2  6z2 


2.2.  Boundary  conditions 

(1)  At  the  bottom  surface  of  hot  plate. 


(3) 


“Cj-  h(Tgas  7) 


(4) 


Table  2 

The  equivalent  thermal— electrical  properties  when  the  TEG  module  with  n  pairs  of 
p-type  and  n-type  semiconductors  is  simplified  to  one  pair. 


Thermo-electrical 

p-Type  and  n-type 

When  TEG  with  n  pairs 

property 

(for  one  pair) 

is  simplified  to  1  pair 

Thermal  conductivity 

kp-N 

kp-N  (Unchanged) 

(W/m  I<) 

Electrical  resistivity 

Pp-n 

k2Pp-n 

(Q  m) 

Seebeck  coefficient 

ap-N 

nap-N 

of  (V/K) 

A  1  ®T 

and  -  kP— 

1  av 


7  ay 

1  d^-J(aP~acond)T  (forp- 


-type) 


and 


dT 


(9) 


(10) 


The  last  terms  in  the  right  side  of  eqs.  (9)  and  (10)  are  the  Peltier 
heat  due  to  the  current  flow  through  two  different  materials  with 
different  Seebeck  coefficients. 


(6)  At  the  top  surface  of  the  TEG. 


-kcer-Q^  —  hc(T  —  Twater)  (11) 

where  hc  and  Twater  represent  the  heat  transfer  coefficient  and 
cooling  water  temperature,  respectively.  Additionally,  the  hc  can  be 
estimated  from  the  Dittus-Boelter  equation,  eq.  (5b).  For  water 
velocity  v  =  0.2  m/s  and  Dh  =  10  mm,  the  hc  is  around  1000  W/m2  K. 
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(7)  At  the  symmetrical  planes  normal  gradients  are  equal  to 
zero.  In  addition,  all  remaining  surfaces  are  set  to  be  adiabatic 
conditions. 

(8)  The  boundary  condition  for  the  voltage  is  V  =  0  V  at  p-type 
thermoelectric  element. 

2.2.1.  Performance  of  a  TEG 

The  performance  of  a  thermoelectric  generator  can  be  charac¬ 
terized.  The  power  generated  from  a  thermoelectric  generator  can 
be  calculated  as  follows: 

Pm  I  V  (12) 

The  TEG  conversion  efficiency  is  given  by 

,  -  £  (.3, 

where  Qh  is  the  heat  transfer  rate  of  the  hot  side 

Oh  =  ~khpA§  (14) 


2.3.  Numerical  methods 

Since  there  is  no  analytical  solution  in  this  study,  the  governing 
equations  eqs.  (1)  and  (3)  along  with  the  boundary  conditions  eqs. 
(4)— (11)  are  solved  numerically  using  the  finite  difference  method. 
Fig.  2  shows  the  physical  model  and  computation  grid  system  for 
the  thermoelectric  module  with  the  spreader  in  Cartesian  coordi¬ 
nates.  A  grid  system  of  83,722  grid  points  was  adopted  typically  in 
the  computation  domain.  Here,  the  grid  pattern  is  coarsely  dis¬ 
played  for  the  convenience  of  the  reader  to  visualize  the  compu¬ 
tational  regions. 

Prior  to  computation,  a  thorough  verification  of  the  grid- 
independence  of  the  numerical  solution  is  performed  in  order  to 
ensure  the  accuracy  and  validity  of  the  numerical  results.  A  grid 
system  of  83,722  grid  points  is  typically  adopted  in  the  computa¬ 
tional  domain.  In  order  to  validate  the  independency  of  the  solution 
on  the  grid,  three  different  grid  systems  are  investigated,  which 
include  about  61,636,  83,722  and  100,954  nodes,  respectively,  for 
a  thermoelectric  module  with  a  spreader.  It  is  found  that  the  relative 
errors  in  the  local  temperature  (T)  and  electric  potential  (V)  between 
the  solutions  of  83,722  and  100,954  are  less  than  3.0%.  Computa¬ 
tions  were  performed  on  an  INTEL  Core  i7-2600  3.40G  personal 
computer,  and  typical  CPU  times  were  about  1500  s  for  each  case. 


Tm„er=  300  K  hc=  1 000  W/m2K 


Ceramic 

substrate 


N 

P 

Conductive 

copper 

r=0V _ 


Spreader 


Teas  It  /. 


physical  model 


b 


computational  grid  system 


Fig.  2.  The  physical  model  and  computational  grid  system. 
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3.  Optimization 

In  the  present  study,  the  simplified  conjugate-gradient  method 
(SCGM)  [26]  has  been  combined  with  finite  differential  method 
(FDM)  code  as  an  optimizer  to  search  the  optimum  TEG  module 
spacing  (S)  and  corresponding  spreader  thickness  (Hsp).  The 
objective  functions  J0bj(x  1**2)  are  defined  as  the  maximum  power 
per  unit  surface  area  (i.e.  the  maximum  power  density,  PmaxIA). 

Above  all,  the  SCGM  method  evaluates  the  gradient  of  the 
objective  function,  and  then  it  sets  up  a  new  conjugate  direction 
for  the  updated  design  variables  with  the  help  of  a  direct  numer¬ 
ical  sensitivity  analysis.  The  initial  guess  for  the  value  of 
each  search  variable  is  made,  and  in  the  successive  steps,  the 


conjugate-gradient  coefficients  and  the  search  directions  are 
evaluated  to  estimate  the  new  search  variables.  The  solutions  ob¬ 
tained  from  the  finite  difference  method  are  then  used  to  calculate 
the  value  of  the  objective  function,  which  is  further  transmitted 
back  to  the  optimizer  for  the  purpose  of  calculating  the  consecu¬ 
tive  searching  directions.  The  procedure  for  applying  this  method 
is  described  in  the  following: 

(1)  Generate  an  initial  guess  for  two  design  variables  (xi,X2)  -  TEG 
module  spacing  (S)  and  its  spreader  thickness  (Hsp). 

(2)  Adopt  the  finite  difference  method  to  predict  the  temperature 
fields  (T)  and  electric  potential  (V)  associated  with  the  latest  S 
and  Hsp,  and  then  calculate  the  objective  function  J0bj(x i,X2). 


Fig.  3.  Flow  charts  for  the  optimization  method. 
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(3)  When  the  value  of  J0bj(x i,x2)  reaches  a  maximum,  the  optimi¬ 
zation  process  is  terminated.  Otherwise,  proceed  to  step  (4). 

(4)  Determine  the  gradient  functions,  {dj0bjldx i)(k)  and  {dj0bjldx 2fk\ 
by  applying  a  small  perturbation  (Axi,Ax2)  to  each  value  of  xi 
and  x2,  and  calculate  the  corresponding  change  in  objective 
function  (A f0bj)-  Then,  the  gradient  function  with  respect  to 
each  value  of  the  design  variables  (xi,x2)  can  be  calculated  by 
the  direct  numerical  differentiation  as 


dJobj  _  A Jobj  djpbj  _  A lobj 

dXj  ~  AX!  0X2  _  Ax2 


(15) 


(5)  Calculate  the  conjugate-gradient  coefficients,  y^  and  and 
the  search  directions,  and  %2+^’  f°r  each  search  variable. 
For  the  first  step  with  k  =  1,  y^  =  0  and  y^  =  0. 


(tylobj/dX  1 

)*i 

-(dJobj/dX  1  ^ 

k-t 

and  y\ 


(pJobj/^2 

)*i 

-(dJobj/dX2^ 

k- 1 

(16) 


axi 


■  Yi? 


and 


pk+ 1 


<bj 

0*2 


+  t& 


(17) 


(6)  Assign  values  to  the  coefficients  of  descent  direction  (/3i,/32)  for 
all  values  of  the  design  variables  (xi,x2).  Specifically,  those 
values  are  chosen  by  a  trial-and-error  process.  In  general,  the 
coefficients  of  descent  direction  (/?i,/?2)  are  within  a  range  of 
0.01-0.001. 

(7)  Update  the  design  variables  with 


xf+1)  =  xf>  +|S1?f)  and  x(,fc+1)  =  x®  4-  M2]  (18) 

A  flowchart  of  the  SCGM  optimization  process  is  plotted  in  Fig.  3. 


module  and  an  exhaust  pipe.  The  air  flow  was  driven  by  a  3.7  kW 
centrifugal  fan  with  an  inverter  to  provide  various  inlet  velocities. 
During  the  isothermal  test,  it  was  determined  that  the  variations  in 
the  temperature  in  different  locations  were  within  ±0.2  °C  and  that 
a  free  stream  turbulence  intensity  of  less  than  2%  was  achieved.  The 
air  temperatures  at  the  inlet  and  the  exit  zones  across  the  test 
section  were  measured  by  two  pre-calibrated  RTDs  (pt-100  Q). 
Their  accuracy  was  within  ±0.05  °C.  The  desired  air  velocity 
through  the  test  section  was  set  by  means  of  adjusting  an  inverter. 
The  air  velocity  was  measured  by  a  hot  wire  with  ±2.0%  accuracy. 
To  reduce  heat  loss,  all  wind  tunnel  system  was  surrounded  by 
insulators. 

The  TEG  module  was  cooled  by  water,  and  its  temperature  was 
controlled  at  27  °C  by  a  thermostat  reservoir.  The  water  volumetric 
flow  rate  was  measured  by  a  magnetic  volume  flow  meter  with 
0.002  L/s  resolution.  The  heat  transfer  coefficients  in  the  water 
channel  (hc)  and  wind  tunnel  (h/0  could  be  estimated  by  using  the 
Dittus-Boelter  equation.  A  high  power  electrical  resistor  array  was 
connected  in  series  to  the  system  which  was  used  to  capture  the 
matching  load  for  the  optimized  power  output.  The  ammeter  and 
the  voltmeter  were  connected  to  measure  the  current  in  the  circuit 
and  voltage  cross  the  external  load  resistor  respectively.  A  digital 
multi-meter  (Fluke  170)  was  used  to  measure  the  voltage  with  an 
accuracy  of  0.15%  and  the  current  with  an  accuracy  of  1.0%.  All  the 
data  signals  were  collected  and  converted  using  a  data  acquisition 
system  (a  hybrid  recorder). 

As  estimating  the  uncertainty  in  measured  and  calculated 
quantities,  the  bias  and  precision  errors  were  both  considered. 
These  elemental  errors  are  combined  to  obtain  an  overall  uncer¬ 
tainty  in  a  measured  quantity  using  the  root-sum-squares  (RSS) 
method.  Uncertainties  in  the  reported  experimental  values  of 
effective  heat  transfer  coefficient  ( h )  and  power  output  (P)  as 
estimated  by  the  method  suggested  by  Moffat  [28].  The  uncer¬ 
tainties  are  tabulated  in  Table  3.  To  minimize  the  precision  uncer¬ 
tainty  of  a  measurement,  the  experiment  was  measured  ten  times 
in  the  same  way  and  in  the  same  circumstance  and  taking  the 
average. 


4.  Experiment  setup  and  uncertainty  analysis 

As  illustrated  in  Fig.  4,  the  wind  tunnel  system  consists  of 
a  blower,  12  annular  heaters  (35  kW),  a  settling  chamber,  the  TEG 


5.  Results  and  discussion 

Fig.  5(a)  and  (b)  illustrates  the  temperature  and  electric  poten¬ 
tial  distributions  of  the  thermoelectric  module  with  and  without 
the  spreader,  respectively,  for  a  TEG  module  spacing  5  =  200  mm, 
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■  Circulating  water 
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Fig.  4.  Schematic  diagram  of  the  experimental  setup. 
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Table  3 

Summary  of  estimated  uncertainties. 


Primary  measurements 

Derived  quantities 

Parameter 

Uncertainties 

Parameter 

Uncertainties 

Uncertainties 

Re  gas  =  1000 

Re  gas  =  5000 

V/r.  gas 

0.3-1% 

Re  gas 

±1.2% 

±0.41% 

Wwater 

0.5% 

Qgas 

±3.5% 

±1.6% 

I 

0.2% 

K 

±5.6% 

±2.2% 

V 

0.1% 

P 

±0.22% 

±0.22% 

Tgas 

0.1  °C 

T water 

0.05  °C 

and  its  spreader  thickness  Hsp  =  5  mm,  where  the  hot  and  cold 
surfaces  were  assumed  to  be  exposed  to  heat  transfer  with 
convective  heat  transfer  coefficient  hh  =  20  W/m2  K  for  the  hot  side 
and  hc  =  1000  W/m2  K  for  the  cold  side  and  with  a  temperature 

a  without  spreader 


difference  (AT  =  200  K).  It  is  seen  from  Fig.  5(a)  that  a  proper  size  of 
heat  spreader  can  indeed  decrease  the  total  thermal  resistance  due 
to  a  larger  heat  conduction  cross  section  area  from  the  hot  gas  side 
to  the  cold  water  side.  Thus,  the  temperature  difference  between 
the  top  and  bottom  surfaces  of  the  TEG  with  a  spreader  is  103  K, 
which  is  larger  than  that  without  a  spreader  (87  K).  This  results  in 
an  open  circuit  voltage  Voc  =  1.133  V  (21%  increase)  for  a  TEG 
module  with  a  spreader  compared  to  Voc  =  0.935  V  for  the  one 
without  a  spreader,  as  seen  in  Fig.  5(b). 

Fig.  6(a)  and  (b)  shows  the  voltage— current  (V-/)  and  power- 
current  (P— J)  curves,  respectively,  for  three  different  values  of 
temperature  difference  between  the  waste  gas  and  the  cooling 
water  (AT  =  100  K,  150  K  and  200  I<)  with  S  =  200  mm,  Hsp  =  5  mm 
and  hh  =  20  W/m2  K.  The  dashed  lines  denote  the  results  without 
a  spreader.  In  addition,  the  experimental  results  are  denoted  by 
circle  symbols.  It  can  be  seen  that  the  predicted  numerical  data  for 
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Fig.  5.  The  temperature  and  electric  potential  distributions  for  the  thermoelectric  module  with  and  without  the  spreader,  respectively,  for  S  =  200  mm,  Hsp  =  5  mm,  hh  =  20  W/m2  K 
and  AT  -  200  K. 


a  1.8 


the  V—I  and  P—I  curves  are  in  good  agreement  (within  8%)  with  the 
experimental  data.  It  can  be  found  that  the  maximum  power,  Pmax, 
occurs  when  the  external  load  resistance  is  to  equal  the  effective 
internal  resistance  of  the  TEG.  A  close  look  at  Fig.  6(b)  indicates  that 
the  Pmax  with  a  spreader  can  be  significantly  increased  up  to  50% 
compared  to  those  without  a  spreader. 

Fig.  7(a)-(c)  presents  the  variations  of  PmaXl  r\  and  Pmax/A  with 
the  TEG  module  spacing  module  spacing  (S),  respectively,  for  four 
different  values  of  waste  gas  heat  transfer  coefficients  (hh  =  20, 40, 


Fig.  6.  V-I  and  P—I  curves  for  different  values  of  temperature  differences  between  the 
waste  gas  and  the  cooling  water  with  S  =  200  mm,  Hsp  =  5  mm  and  hh  =  20  W/m2  K. 


power  output  ( Pmax ) 


S  (mm) 

power  density  (Pma*IA) 

Fig.  7.  The  Pmax,  7]  and  Pmax/A  vs.  TEG  module  spacing  (S)  for  different  values  of  waste 
gas  heat  transfer  coefficients  with  Hsp  =  5  mm  and  AT  =  200  K. 
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S  (mm) 

power  density  (PmskX/A) 


60  and  80  W/m2  I<)  with  a  spreader  thickness  Hsp  =  5  mm  and 
AT  =  200  K.  As  expected,  both  the  power  output  and  TEG  efficiency 
(rj)  are  increased  as  either  the  module  spacing  (S)  or  the  waste  gas 
heat  transfer  coefficient  (hh)  is  increased  (this  corresponds  to 
a  larger  heat  transfer  area  A  =  S2  or  a  higher  waste  gas  velocity). 
However,  there  is  an  asymptotic  nature  of  the  Pmax-S  curves  when 
S  is  greater  than  around  150  mm.  Thus,  there  is  optimum  module 
spacing  S  for  the  maximum  power  density  (Pmox/A),  as  seen  from 
Fig.  7(c).  As  an  example,  for  Hsp  =  5  mm  and  AT  =  200  K,  the  optimal 
module  spacing  S  is  about  110  mm  for  hh  =  20  W/m2  K,  while  the 
optimum  S  is  about  60  mm  for  hh  =  80  W/m2  K. 

Fig.  8(a)— (c)  demonstrates  the  variations  in  Pmax,  rj  and  Pmax/A 
with  the  TEG  module  spacing  (S)  for  four  temperature  differences 


Fig.  8.  The  Pmax,  7]  and  Pmax/y\  vs.  TEG  module  spacing  (S)  for  different  values  of 
temperature  differences  between  the  waste  gas  and  the  cooling  water  with 
Hsp  =  5  mm  and  hh  =  20  W/m2  K. 


Fig.  9.  The  Pmax,  PmaxlA  and  tj  vs.  spreader  thickness  (Hsp)  for  different  values  of  waste 
gas  heat  transfer  coefficients  with  S  =  200  mm  and  AT  =  200  K. 
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between  the  waste  gas  and  the  cooling  water  (AT  =  200,  400,  600 
and  800  I<)  with  spreader  thicknesses  of  Hsp  =  5  mm  and 
h/i  =  20  W/m2  K.  It  is  evident  that  both  the  power  output  and  TEG 
efficiency  (rj)  are  significantly  increased  as  the  temperature 
difference  AT  between  the  waste  gas  and  the  cooling  water  is 
increased.  However,  the  optimum  module  spacing  (5)  for  the 
maximum  power  density  ( PmaxIA )  is  weakly  dependent  on  AT.  For 
example,  for  Hsp  =  5  mm  and  hh  =  20  W/m2  K,  the  optimal  module 
spacing  S  is  about  116  mm  for  AT  =  200  K,  while  the  optimum 
S  =  112  mm  for  AT=  800  K. 

Figs.  9  and  10  show  the  effect  of  the  spreader  thickness  (Hsp)  on 
Pmax  and  I)  for  different  values  of  hh  and  AT,  respectively,  with 
S  =  200  mm.  It  can  be  seen  that  the  power  is  first  increased  as  the 


power  output  (Pmax)  and  power  density  (Pm.JA) 
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Fig.  10.  The  Pmax,  PmaxIA  and  tj  vs.  spreader  thickness  (Hsp)  for  different  values  of 
temperature  differences  between  the  waste  gas  and  the  cooling  water  with 
S  =  200  mm  and  hh  =  20  W/m2  K. 
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Fig.  11.  Iteration  process  to  search  the  optimum  combination  of  S  and  Hsp  for  the 
maximization  of  objective  function  with  hh  =  20  W/m2  K  and  AT  =  200  K. 


spreader  thickness  {Hsp)  is  increased.  This  is  due  to  the  fact  that,  as 
the  spreader  thickness  is  increased  at  the  initial  stage,  there  is  an 
advantage  of  descending  spreading  resistance  due  to  the  path 
associated  with  the  lateral  heat  diffusion  from  the  hot  waste  gas  to 
the  TEG  module.  However,  as  the  spreader  thickness  (Hsp)  is 
increased  over  about  60  mm,  there  is  greater  thermal  resistance 
penalty  over  the  reduced  spreading  resistance,  and  thus,  the  power 
may  be  decreased.  The  TEG  conversion  efficiency  (rj)  is  increased  as 
spreader  thickness  (Hsp)  is  increased.  However,  there  is  an 
asymptotic  nature  of  the  r]-Hsp  curves  when  Hsp  is  greater  than 
around  10  mm. 

Figs.  11-13  display  three  examples  of  the  iteration  process  to 
search  the  optimum  TEG  module  spacing  (S)  and  its  spreader 
thickness  (Hsp)  combination  for  the  maximization  of  objective 
function  (i.e.  power  density,  PmaxlA)  for  three  operational 


Fig.  12.  Iteration  process  to  search  the  optimum  combination  of  S  and  Hsp  for  the 
maximization  of  objective  function  with  hh  =  60  W/m2  K  and  AT  =  200  K. 
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Fig.  13.  Iteration  process  to  search  the  optimum  combination  of  S  and  Hsp  for  the 
maximization  of  objective  function  with  =  20  W/m2  K  and  AT  =  600  K. 


conditions  (1)  hh  =  20  W/m2  K,  AT  =  200  K  (2)  hh  =  60  W/m2  K, 
A T  =  200  K  (3)  hh  =  20  W/m2  K,  AT  =  600  K,  respectively.  The 
constant  power  density  contours  are  plotted  as  a  function  of  S  and 
Hsp ,  where  the  dark  brown  area  represents  the  maximum  power 
density.  It  is  noted  that  these  constant  power  density  contours  were 
obtained  by  using  400  computer  simulations  for  each  plot.  It  is  seen 
that,  by  using  the  simplified  conjugated-gradient  method  (SCGM), 
the  optimal  S  and  Hsp  combination  can  be  obtained  only  for  around 
40  iterations  (case  1:  31  iterations,  case  2:  35  iterations,  case  3:  41 
iterations).  Thus,  the  current  optimization  method  provides 
a  tremendous  savings  to  computational  time  in  the  present  study. 

The  searched  optimum  combination  of  S  and  Hsp  for  different 
operational  conditions  (AT  =  200  -  800  K  and  hh  =  20  -  80  W/ 
m2  K)  is  tabulated  in  Table  4.  It  can  be  seen  that  the  optimum  (S,  Hsp) 
are  strongly  dependent  on  the  waste  gas  heat  transfer  coefficient, 
hh ,  and  they  are  decreased  as  hh  is  increased,  while  they  are  weakly 
dependent  on  the  temperature  difference  between  the  waste  gas 
and  the  cooling  water,  AT. 


Table  4 

The  searched  optimum  combination  of  S  and  Hsp  for  different  operational  conditions. 


AT(K) 

hh  (W/m2  K) 

S  (mm) 

Hsp  (mm) 

200 

20 

116.5 

21.6 

40 

82.3 

13.3 

60 

66.8 

9.4 

80 

58.1 

7.1 

400 

20 

114.6 

20.9 

40 

80.9 

13.1 

60 

66.0 

9.3 

80 

57.4 

7.0 

600 

20 

113.2 

20.6 

40 

79.9 

12.9 

60 

65.3 

9.1 

80 

56.8 

6.9 

800 

20 

112.4 

20.3 

40 

79.4 

12.8 

60 

64.8 

9.0 

80 

56.5 

6.8 

6.  Conclusions 

In  this  study,  first,  the  parametric  study  for  the  thermal-electrical 
characteristics  of  a  TEG  module  used  in  a  waste  heat  recovery  system 
was  investigated  numerically  and  experimentally.  It  was  demon¬ 
strated  that  the  proper  size  of  a  heat  spreader  can  decrease  the 
thermal  resistance  and  that  the  maximum  power  Pmax  with 
a  spreader  can  be  significantly  increased  (up  to  50%)  as  compared  to 
TEG  without  spreader.  The  predicted  numerical  data  for  the  power  vs. 
current  (P— /)  curve  are  in  good  agreement  (within  8%)  with  the 
experimental  data.  Secondly,  the  optimization  of  the  TEG  module 
spacing  (S)  and  corresponding  spreader  thickness  (Hsp)  was  per¬ 
formed  by  using  a  simplified  conjugate-gradient  method.  The  power 
density  was  the  objective  function  to  be  maximized.  The  searched 
optimum  combination  of  S  and  Hsp  for  different  operational  condi¬ 
tions  (AT  =  200  -  800 1<  and  hh  =  20  -  80  W/m2 1<)  was  obtained.  The 
optimum  TEG  module  spacing  (S)  and  its  spreader  thickness  (Hsp) 
were  found  to  be  strongly  dependent  on  the  waste  gas  heat  transfer 
coefficient,  h h\  while  they  were  weakly  dependent  on  the  tempera¬ 
ture  difference  between  the  waste  gas  and  the  cooling  water,  AT.  The 
results  obtained  in  this  study  should  provide  useful  information  for 
the  industrial  TEG  module  design  for  waste  heat  recovery. 
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